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Design of a Mott Multiferroic from a Non-Magnetic Polar Metal 
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We examine the electronic properties of newly discovered “ferroelectric” metal Li0s03 combining 
density-functional and dynamical mean-field theories. We show that the material is close to a Mott 
transition and that electronic correlations can be tuned to engineer a Mott multiferroic state in 1/1 
super lattice of LiOsOs and LiNbCW We use electronic structure calculations to predict that the 
(Li0s03)i/(LiNb03)i superlattice is a type-I multiferroic material with a ferrolectric polarization of 
41.2 pC cm -2 , Curie temperature of 927K, and Neel temperature of 671 K. Our results support a 
route towards high-temperature multiferroics, z.e., driving non-magnetic polar metals into correlated 
insulating magnetic states. 
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Introduction .—Multiferroics (MF) are a class of insu¬ 
lating materials where two (or more) primary ferroic 
order parameters, such as a ferroelectric polarization 
and long-range magnetic order, coexist. Technologically, 
they offer the possibility to control magnetic polariza¬ 
tions with an electric field for reduced power consumption 
[1, 2]. Nonetheless, intrinsic room-temperature MF re¬ 
main largely elusive. This fact may be understood by ex¬ 
amining the microscopic origins for the ferroic order which 
aids in classifying different phases: In Type-I MF, ferro- 
electricity and magnetism arise from different chemical 
species with ordering temperatures largely independent 
of one another and weak magnetoelectric (ME) coupling 
[3]. The ferroelectric ordering also typically appears at 
temperatures higher than the magnetic order, and the 
spontaneous polarization P is large since it is driven by 
a second-order Jahn-Teller distortion, e.g ., BiFeOs [3, 4]. 
In Type-II MF, however, magnetic order induces ferro- 
electricity, which indicates a strong ME coupling between 
the two order parameters. Nonetheless, P is usually much 
smaller, e.g., by a factor of 10 2 as in R-M^Os (R be¬ 
ing rare earth) [5]. In a few MFs with high-transition 
temperatures, z.e., BiFeOs [6] and Sri-^Ba^MnOs [7-9], 
magnetism is caused by Mott physics arising from strong 
correlations. The interactions localize the spins at high 
temperature, paving the way for magnetic ordering at 
room temperature. Materials where this robust mag¬ 
netism is coupled with ferroelectric distortions are ideal 
candidates for a room-temperature MFs. 

Herein, we propose a design strategy for novel Mott 
MF phases. It relies on tuning the degree of correlation 
of the recently discovered class of materials referred to 
as ‘ferroelectric metals’ with LiOsOs as the prototypi¬ 
cal member [10]. This material is the first undisputed 
realization of the Anderson-Blount mechanism [11], and 
challenges the expectation that conduction electrons in 
metals would screen the electric field induced by polar 
displacements [10, 12, 13]. Despite robust metallicity, this 
material shares structural similarities with prototypical 
insulating ferroelectric oxides, such as LiNbOs [14, 15]: A 


R3c crystal structure with acentric cation displacements 
and distorted 0s06 octahedra [16, 17] and comparable 
lattice parameters [10, 14]. While the polar displacements 
in LiNbOs rely on cross-gap hybridization between p (O) 
and d (Nb) states [18], in LiOsOs they are weakly coupled 
to the states at the Fermi level (£/?), which makes possi¬ 
ble the coexistence of an acentric structure and metallicity 
[16, 19]. In LiOsOs the empty d-manifold of LiNbOs is 
replaced by a non-magnetic 5 d 3 ground state with a half- 
filled t 2 g ( d xy , d xz , d yz ) configuration, which is responsible 
for the metallic response [16]. However, the strength of 
the electronic interactions is insufficient to drive a Mott 
transition in the correlated t^g manifold as revealed by 
low-temperature resistivity measurements; nonetheless, if 
it would be possible to enhance the electronic correlations 
in LiOsOs and achieve a metal-insulator transition, then a 
previously unidentified multiferroic material should result. 
The concept is that if an insulating state can be obtained 
from a ‘ferroelectric metal’ through enhanced correlations, 
it would then naturally lead to magnetic ordering of the 
localized electron spins, coexisting polar displacements, 
and potentially strong ME coupling. 

In this work we explore the feasibility of this ap¬ 
proach using a combination of first-principles density 
functional theory (DFT) plus dynamical mean field the¬ 
ory (DMFT) calculations [20]. We first show that the 
electronic Coulomb interactions and Hund’s coupling in 
LiOsOs make it an ideal candidate for realizing a Mott 
MF due to the multi-orbital t^g physics. Next, we de¬ 
scribe the design of a new multiferroic by control of the 
electronic structure through atomic scale engineering of a 
Mott metal-insulator transition (MIT) in an ultrashort 
period (Li 0 s 03 )i/(LiNb 03 )i superlattice. The insulat¬ 
ing and magnetic state is driven by an enhancement of 
the electronic correlations in LiOsOs layers owing to the 
kinetic energy reduction of the t 2 g orbitals from the su¬ 
perlattice geometry. The ferroelectric properties mainly 
originate from cooperative Li and O displacements. The 
multiferroic phase emerges across the MIT, exhibiting 
a net electric polarization (41.2 pC cm -2 ) and magneti- 
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zation [0.9/i# per formula unit (f.u.)], with calculated 
magnetic-ordering and ferroelectric temperatures of 671 K 
and 927 K, respectively. Our results uncover a promising 
alternative route to discovery of room-temperature mul- 
tiferroics: One could search for correlated polar metals 
near Mott transitions and drive the phases into insulating 
states, rather than the often-pursed approach of inducing 
polar displacements in robustly insulating magnets. 

Calculation Details. —We perform first-principles DFT 
calculations within local-density approximation (LDA) 
+Hubbard U method as implemented in the Vienna A b 
initio Simulation Package (VASP) [21] with the projector 
augmented wave (PAW) approach and a 600 eV plane 
wave cutoff with a 5 x 7 x 7 k- point mesh. We relax the 
volume and atomic positions (forces <0.1 meV A -1 ) using 
Gaussian smearing (20meV width) for the Brillouin zone 
(BZ) integrations. We perform LDA+DMFT calculations 
including local Coulomb interactions parameterized by 
the U and the Hund’s coupling J\ starting from Wannier 
orbitals constructed from the LDA bands [22] using an 
energy range spanned by the full d manifold. The impurity 
model is solved using Exact Diagonalization (ED) with a 
parallel Arnoldi algorithm [23, 24]. 

Correlations in LiOsOs. —We first examine the effect 
of the interactions on the metallic state of LiOsOs and 
determine the critical values for a Mott transition U c in 
the paramagnetic and antiferromagnetic (AFM) phases 
using LDA+DMFT. The criterion for a Mott-Hubbard 
transition is frequently associated with the ratio between 
the bandwidth (W) and the interaction strength V, so 
that the Mott transition occurs for U c of the order of W. 
In a multiband Hubbard model with M orbitals, U c is 
enhanced by orbital fluctuations, i.e., U c VMW, [25] 
and it is influenced by the Hund’s coupling J^. Indeed, at 
half-filling, U c is reduced by an enhancement of Jh [26]. 

In the following, we show this is precisely the situation 
in LiOsOs [16]. Due to the energy separation between t^g 
and e g orbitals in the spectral density of state of LiOsOs 
around the Fermi level, we resort to using a model for 
the t 2 g levels only [16]. Owing to the symmetry breaking 
in bulk LiOsOs, the orbitals in the d manifold are also 
permitted to mix, which lifts the degeneracy of t^ g orbitals 
with two of states remaining degenerate. 

Fig. 1 shows the orbital resolved quasiparticle weight 
(Z) of the occupied orbitals as a function of U for two 
different values of Jh for paramagnetic LiOsOs in the 
experimental structure (see top panels). Z measures the 
metallic character of the system, and it evolves from Z=1 
for a non-interacting metal to Z=0 for a Mott insulator. 
Upon increasing the value of J^, the critical value of U 
required to reach the Mott state (Z=0) is shifted to larger 
values of U [26] . 

In the correlated regime, we anticipate electron localiza¬ 
tion will lead to long-range magnetic order of the localized 
spins. Spin-polarized LDA+DMFT calculations, initial¬ 
ized with a G-type AFM structure (every spin on an Os 
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FIG. 1. (Color online) (a) Orbital resolved quasiparticle orbital 
weight Z (filled symbols) for paramagnetic LiOsOs and (b) 
local magnetization m (/jlb) (obtained from a spin-polarized 
calculation) of the t 2 g orbitals in LiOsOs as function of U for 
different ratios of Jh/U within the LDA+DMFT calculations. 
Vertical arrows indicate the critical value of U required to 
reach the insulating state in the G-type AFM structure. 

cation is antiparallel to all its neighbors), reveal that the 
local magnetic moment rapidly saturates to the atomic 
value S = 3/2. A finite magnetization also develops at 
intermediate U in the metallic state (Fig. 1, lower panels). 
The MIT, marked by vertical arrows, occurs for a weaker 
coupling in the AFM than in the paramagnetic state. 

Design of a Mott Multiferroic. —The LDA+DMFT cal¬ 
culations reveal that a simultaneous Mott and magnetic 
state could be engineered in LiOsOs by reducing the 
electronic kinetic energy. One avenue to control and 
decrease the kinetic energy relies on heterostructuring 
and interleaving two perovskites together to form a co¬ 
herent superlattice, whereby an isostructural insulator 
would restrict the electron hopping due to the reduction 
in available channels [27-29]. Such geometries can be 
achieved in practice using oxide molecular-beam epitaxy 
or pulsed-laser deposition methods [30, 31]. 

Owing to the chemical and structural compatibility of 
LiOsOs with LiNbOs, with a lattice mismatch of 3.2%, 
we devise an ultrashort period perovskite superlattice 
of (LiOsOs)i/(LiNbOs)i as illustrated in panel (a) of 
Fig. 2. The superlattice is constructed by beginning from 
the R3c crystal structure of LiOsOs (LiNbOs) and im¬ 
posing a layered order along the [110] direction in the 
rhombohedral setting, which is equivalent to a 1/1 period 
LiOsOs/LiNbOs grown along the pseudocubic (pc) [001] 
direction [32]. The geometry in Fig. 2 is also different from 
a superlattice constructed along the [101] pc direction (c./., 
Ref. 33), which is likely more challenging to realize ex¬ 
perimentally. Following full relaxation of the superlattice, 
we find the cation order results in a symmetry reduction 
to the polar space group Pc with out-of-phase Os06 and 
NbOe octahedral rotations, z.e., the a~b~b~ tilt pattern 
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FIG. 2. (Color online) (a) The superlattice exhibits the 
a~b~b~ tilt pattern. Atom- and orbital-resolved DOS for 
(b) Li 0 s 03 , (c) LiNb 03 and (d) LiOsOs/LiNbOs at the DFT- 
LDA level. Ef is given by the (broken) vertical line at OeV. 


given in Glazer notation [34]. The microscopic origin of 
polar displacements are described in detail below. 

Electronic Properties. — Fig. 2 shows the LDA electronic 
density of states (DOS) for the LiOsOs/LiNbOs superlat¬ 
tice (d), compared with LiOsOs (b) and LiNbOs (c) using 
the LDA-optimized atomic structures. The results for 
LiOsOs (Fig. 2b) highlight the metallic character of the 
former, where the weight at the Fermi level (Ef) mainly 
comes from Os 5 d states which show strong admixture 
from the O 2 p states. In contrast, LiNbOs is a band 
insulator, with the O 2 p states forming the valence band 
and Nb 4 d states at the conduction band minimum, sepa¬ 
rated by a gap of 3.28 eV (Fig. 2c). In the superlattice, 
we find essentially no charge transfer between Os and Nb: 
Each component (LiOsOs and LiNbOs) is isoelectronic 
to its bulk configuration; the DOS can be described as a 
direct superposition of the two components (Fig. 2d). The 
Os 5 d states partially fill the gap in the electronic spec¬ 
trum formed from the the two-dimensional Nb 02 planes. 
There is some spectral weight transfer in the vicinity of 
Ep among the Os orbitals, which are sensitive to the 
electron correlation strength as shown in Fig. 1. 

We now explore the effect of electronic correlations 
by means of LSDA+L/ calculations at different values of 
£4ff = U — Jh- An accurate value of the Hubbard U is 



FIG. 3. (Color online) Band gap E g and averaged local mag¬ 
netic moment for Os as a function of C7 e ff with and without 
spin-orbit interaction (SOI). 


unknown for perovskite osmates, but it is expected to be 
comparable to that of NaOsOs [35] and double perovskite 
Sr2Cr0s06 [36] for which a correct description of the 
electronic properties are obtained with U values of 1.0 
and 2.0 eV, respectively. Note that the differences from 
various implementations of the LDA -\-U scheme for bulk 
LiOsOs were found to be minor [16], and are anticipated 
to also be insignificant for the superlattice. 

Fig. 3 shows the evolution in the band gap (E g ) and 
magnetic moment of Os 3+ ions (m) as a function of the 
strength of U e s for LSDA including spin-orbit interaction 
(SOI, broken lines). A gap opens at a critical [/ e ff ~ leV 
(£/c), signaling a MIT into a magnetic insulating ground 
state. As expected the LiOsOs/LiNbOs superlattice be¬ 
comes insulating for smaller values of the interaction with 
respect to bulk LiOsOs. The enhancement of electronic 
correlations is also found for small values of U. In fact, 
LiOsOs remains paramagnetic [16] while the superlattice 
is weakly ferrimagnetic already below L/ c , where it turns 
into a G-type AFM insulator at C7 e ff ~ 0.5 eV. 

The reduction in U c for the MIT in the superlattice 
can be understood by analyzing the effect of the geo¬ 
metrical confinement on the t^g band dispersions. (For 
simplicity, we use the LDA electronic structures given 
in Ref. 32.) While the bandwidth of the d xy orbitals is 
essentially the same as for bulk LiOsOs, the d xz and d yz 
bands in LiOsOs/LiNbOs are significantly narrowed as a 
consequence of the reduced hopping along the superlattice 
direction. This leads to a reduction of the kinetic energy 
which enhances the electron-electron correlations, making 
the superlattice a Mott insulator at moderate interaction 
strengths. We note that when SOI are excluded in the 
calculations (Fig. 3, solid lines), the MIT occurs at a fur¬ 
ther reduced correlation strength (U c ~ 0.5eV), and the 
magnetic moment only slightly increases. Such behaviors 
are also observed in bulk LiOsOs [16]. 

Ferroelectric Polarization .—We now apply a group the- 
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FIG. 4. (Color online) Illustration of the polar zone-center 
mode along the [101]-direction labeled by irrep ■ Anti-polar 
displacements along the [010]-direction are omitted for clarity. 

oretical analysis [37, 38] of the LiOsOs/LiNbOs structure 
to understand the inversion symmetry-breaking displace¬ 
ments that produce the Pc ground state. We use a fic¬ 
titious P2\/c centrosymmetric phase (where polar dis¬ 
placements are switched off) as the reference phase from 
which the symmetry-adapted mode displacements are ob¬ 
tained as different irreducible representations (irreps) of 
the P2i/c space group operators [39]. We find the loss 
of inversion symmetry mainly derives from cooperative 
Li and O displacements in the (101) mirror plane of the 
Pc phase. Moreover, we find anti-polar displacements 
along the 6-axis which result in no net polarization. All 
polar displacements are described by a distortion vector 
that corresponds to the irrep along the [101]-direction 
of the Pc structure (Fig. 4). These displacements are 
consistent with the acentric Li and O ionic displacements 
identified to be responsible for lifting inversion symme¬ 
try in bulk LiOsOs [10, 13] and across the ferroelectric 
transition in LiNbOs [40]. 

We now compute the ferroelectric polarization in 
LiOsOs/LiNbOs using the Berry’s phase approach [41] 
within LSDA+£7 (U c = 0.5 eV). The spontaneous elec¬ 
tric polarization of the Pc phase is 32.3 fiC cm -2 and 
25.5/iC cm -2 along the [100]-direction, i.e., along the 
pseudo-cubic [001] superlattice repeat direction and [001]- 
directions, respectively. (Note that the [101]-direction in 
LiOsOs/LiNbOs corresponds to the polar [11 Indirection 
in LiNbOs.) Together this yields a net polarization along 
the [101]-direction of 41.2 /iCcm -2 . These values are also 
robust to SOI, with a change of less than 15% to value 
of the total polarization. Following the recipe of Ref. 42, 
we use the energy difference between the high-symmetry 
(P2i/c) and low-symmetry (Pc) to obtain a ferroelectric 
Curie temperature of 927 K for the superlattice. This 
value is close to the extrapolated transition temperature 
for LiNbOs (>1,400 K) [43], and far exceeds that of bulk 
LiOsOs where inversion symmetry is lost near 140 K [10]. 

Magnetic Ordering Temperature .—Our DMFT calcula¬ 


tions indicate that when the superlattice material enters 
in the Mott state the magnetic moment is ^3/i#, corre¬ 
sponding to a high-spin S = 3/2 state. We now estimate 
the Neel temperature for LiOsOs/LiNbOs by extracting 
the exchange interaction constants from spin-polarized 
DFT energies computed at U c without SOI following 
the approach in Ref. 44. Assuming that the magnetism 
arises by ordering such localized spins, we obtain intra- 
and inter-plane Os-Os exchange magnetic couplings of 
-5.6 meV and -0.2 meV respectively, where a negative in¬ 
teraction indicate AFM exchange. From these values and 
without Anderson’s renormalization [45], we estimate a 
Neel temperature of 671 K for the LiOsOs/LiNbOs super¬ 
lattice, which makes the material a correlation-induced 
room-temperature multiferroic. 

Conclusions .—We used a LDA+DMFT approach to 
study the electronic properties of the “ferroelectric” metal 
LiOsOs. A detailed understanding of the electronic struc¬ 
ture of LiOsOs shows that a reduction of the kinetic 
energy can drive the system into a Mott insulating state. 
We use this concept to propose a strategy to design mul¬ 
tiferroic materials by constructing a superlattice with 
the uncorrelated polar LiNbOs dielectric. On the basis 
of LSDA+U calculations we show that the ultra-short 
period LiOsOs/LiNbOs superlattice should be a type-I 
room-temperature Mott multiferroic with a large 41.2 /iC 
cm -2 electric polarization. 

The large ferroelectric displacements from the LiNbOs 
layers facilitate the high ferroelectric ordering temperature 
in the LiOsOs/LiNbOsheterostructure as observed from 
the similarity in the Curie temperature of the superlattice 
with that of LiNbOs. I n this case LiOsOs/LiNbOswould 
behave as a paramagnetic Mott ferroelectric at high tem¬ 
peratures and transition into Mott multiferroic below the 
Neel temperature, which is predicted to be well-above 
room temperature. Because the exchange interactions of 
Os are mediated by the coordinating O ligands, which 
are essential to and produce the ferroelectric distortion, a 
strong ME coupling is anticipated as in Sri-^Ba^MnOs 
[7-9]. We hope this work motivates the synthesis of new 
artificial multiferroics, and the adds to the growing dis¬ 
cussion of new applications where noncentrosymmetric 
metals and ferroelectric materials may be united. 
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TABLE I. Calculated crystallographic parameters for R3c LiOsOs using LDA functional. 

LiOsOs 



a — b 

= 4.9399 A, c = 13.3019 A 

R3c 




a = /3 = 90°, 7 = 120° 

Atom 

Wyck. Site 

X 

y 

z 

Li 

6 a 

0 

0 

0.61117 

Os 

6b 

0 

0 

0.32084 

0 

1 8b 

-0.00388 

0.36259 

0.06825 


TABLE II. Calculated crystallographic parameters for R3c LiNb 03 using LDA functional. 

LiNbOs 



a — b — 

5.1052 A, c= 13.7472 A 

R3c 




a = /3 = 90°, 7 = 120 ° 

Atom 

Wyck. Site 

X 

y 

z 

Li 

6a 

0 

0 

0.61517 

Nb 

6b 

0 

0 

0.33169 

O 

18b 

-0.01367 

0.35987 

0.06330 


TABLE III. Calculated crystallographic parameters for Pc L^NbOsOe using the LDA functional. 

LiOsOs/LiNbC>3 

Pc 

Atom 

Wyck. Site 

X 

a 

y 

= 7.35270 A, b = 5.02739 A, c = 5.36541 A 
a = 7 = 90°, /3 = 95.109° 

z 

Lil 


2 a 

0.27896 

0.25751 

0.45909 

Li2 


2 a 

0.78342 

0.74972 

-0.07013 

Nb 


2 a 

-0.00633 

0.25287 

0.00169 

Os 


2 a 

0.49713 

0.75123 

0.53150 

01 


2 a 

0.23965 

0.12376 

0.02510 

02 


2 a 

0.73578 

0.60514 

0.54599 

03 


2 a 

-0.09599 

0.05685 

0.71217 

04 


2 a 

0.41768 

0.54538 

0.23192 

05 


2 a 

0.04578 

0.58866 

0.85929 

06 


2 a 

0.54412 

0.04168 

0.32305 


TABLE IV. Calculated crystallographic parameters for Pc LLNbOsOe without monoclinic angle using the LDA functional. 

LiOsOs/LiNbOs 

Pc 

Atom 

Wyck. Site 

X 

a 

y 

= 7.35884 k,b = 5.03159 A, c = 5.36989 A 
a = 0 = 7 = 90° 

z 

Lil 

2 a 

0.27310 

0.26765 

0.47516 

Li2 

2 a 

0.78081 

0.75357 

-0.05429 

Nb 

2 a 

-0.00421 

0.25110 

0.00188 

Os 

2 a 

0.49815 

0.75234 

0.52637 

01 

2 a 

0.24071 

0.12246 

0.01238 

02 

2 a 

0.73609 

0.60651 

0.54411 

03 

2 a 

-0.09326 

0.05556 

0.71992 

04 

2 a 

0.41773 

0.53696 

0.24015 

05 

2 a 

0.04603 

0.58671 

0.85024 

06 

2 a 

0.54504 

0.03169 

0.30376 



TABLE V. Calculated crystallographic parameters for Pc L^NbOsOe without monoclinic angle using the LSDA+f/ with 
U=U C = 0.5 eV. 


LiOsOs/LiNbOs 

Pc 

Atom 

Wyck. Site 

X 

a 

y 

= 7.35714 A, b = 5.03042 A, c = 5.36865 A 
a = p = 7 = 90° 

z 

Lil 

2 a 

0.27431 

0.26460 

0.47149 

Li2 

2 a 

0.77680 

0.75646 

-0.05271 

Nb 

2 a 

-0.00460 

0.25150 

0.00212 

Os 

2 a 

0.49708 

0.75063 

0.52604 

01 

2 a 

0.24114 

0.12241 

0.01487 

02 

2 a 

0.73596 

0.60719 

0.54175 

03 

2 a 

-0.09280 

0.05536 

0.72002 

04 

2 a 

0.41875 

0.54274 

0.23604 

05 

2 a 

0.04702 

0.58598 

0.84961 

06 

2 a 

0.54653 

0.03942 

0.31046 


CD 

Q_ 

> 

JO) 

^0 

CD 

as 

GO 


</) 

o 

Q 




"I-1-r 


- 2-10 12 3 

Energy (eV) 


FIG. 5. (Color online) Resolved 5 d t 2 g states of Os for (a) LiOsOa and (b) LENbOsOe within LDA. The bandwidth of the 
d xy orbitals for bulk LiOsC >3 and LENbOsOe are the same while in the superlattice the d xz and d yz orbitals have a reduced 
bandwidth. 
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FIG. 6. (Color online) Superlattice construction: (a) Starting from the R3c crystal structure of LiNb 03 in the hexagonal setting 
with 30 atoms, (b) we transform to the rhombohedral setting with 10 atoms, (c) Next we double the cell with the transformation 
matrix and (d) substitute a Nb02 layer with 0s02 layer, (e) We identify the new space group Pc with FINDSYM [1] and finally 
(f) we set the monoclinic angle to 90°. 











